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ABSTRACT 

The results of measurements and computer simula- 
tions of radioactivities and dose rates as functions of 
decay time are presented for not Pb and 209 Bi irradi- 
ated by 1.5-GeV protons, 59 Co, 63 Cu, and 65 Cu irra- 
diated by 0.13- and 1.2-GeV protons, and 232 Th and 
natjj j rrac iiated by 0.1- and 0.8-GeV protons. The ac- 
tivities and dose rates are measured by direct high- 
precision 7-spectrometry. The irradiations were made 
using external beams extracted from the ITEP U-10 
proton synchrotron. Simulations made using the LCS 
and CINDER'90 code systems are compared with mea- 
surements. 

INTRODUCTION 

The design of a hybrid system driven via a high cur- 
rent accelerator needs to have information about resid- 
ual nuclides that are produced by high energy protons 
in the target and in constructional elements, defining a 
number of technical features of the system. 

Our previous works [1,2] present the results of exper- 
imental measurements of residual radioactive nuclides 
produced in the target and constructional materials ir- 
radiated with 0.13- and 1.5-GeV protons. Predictive 
powers of a number of high energy codes were inves- 
tigated using the experimental data. A general con- 
clusion found insufficient accuracy of the codes to sim- 
ulated radionuclide yields, both independent and cu- 
mulative, with the mean deviation of simulated results 
from experimental data often of a factor of 2 or higher. 
Reliable information on the yield of each residual nu- 
clide is very important and such characteristic can be 
defined as "individual". 

Unlike our previous works, this study is aimed to in- 
vestigate the predictive power of codes not for individ- 
ual but for "integral" parameters - time dependent ac- 
tivation and dose-equivalent rate (called dose rate here- 



after) which are directly dependent on summed yield of 
individual radionuclides. For such "integral" character- 
istics, one may expect a much better predictive power 
of codes than in simulations of "individual" parame- 
ters. We believe that the dose rate due to emulations 
of the 7-s from any hybrid system part being irradi- 
ated is one of the most important working parameters 
in applications. 

Here, the LAHET [3] and HMCNP Monte Carlo 
transport codes were used together with the code 
CINDER'90 [4] to simulate the build-up and decay of 
the products. 

EXPERIMENT 

The irradiations of experimental samples were car- 
ried out using external beams of the ITEP proton syn- 
chrotron. A sandwich of experimental and aluminum 
samples of the same diameters of 10.5 mm was placed 
perpendicularly to the proton beam during each irra- 
diation. The parameters of irradiations and of experi- 
mental samples are presented in Table 1. 

The CANBERRA spectrometer based on the GC- 
2518 Ge detector was used to measure 7-spectra of irra- 
diated experimental and aluminum samples. 7-spectra 
processing and nuclide identification were made using 
the ASPRO and SIGMA codes and the PCNUDAT nu- 
clear database. 24 Na formation in aluminum samples 
was used to determine the irradiation time integrated 
proton flux. The details of irradiation techniques, 
schemes and parameters of external beams, techniques 
of nuclides identification and their cross section deter- 
mination can be found in [1]. 

SIMULATION TECHNIQUE 

The simulation of build-up during irradiations and 
following decay of produced nuclides was made using 
the LCS [3] and CINDER'90 [4] code systems. 
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Table 1: Parameters of irradiations and experimental 
samples 



N 


Target 


Parameter 






Sample 


Proton 


Irradi- 


Proton 






weight 


en- 


ation 


flux 






(mg) 


ergy 
(MeV) 


time 
(min) 


(p/cm2) 

during 

irradiation 


A 


natp^ 


143.8 


1478 


15 


2.43E+13 


B 


209 Bi 


368.05 


1478 


15 


2.67E+13 


C 


59 Co 


202.7 


1186 


45 


6.86E+13 


D 


59 Co 


196.6 


127 


60 


2.29E+13 


E 


63 Cu 


80.2 


1186 


70 


1.08E+14 


F 


65 Cu 


92.6 


1186 


70 


1.37E+14 


G 


63 Cu 


425.55 


127 


9 


1.56E+11 


H 


65 Cu 


521.6 


127 


9 


1.49E+11 


I 


232 Th 


87.5 


97 


40 


6.31E+12 


J 


232 Th 


89.6 


795 


30 


4.86E+12 


K 


nat^j 


159.2 


97 


60 


9.22E+12 


L 


natjj 


160.9 


795 


60 


7.43E+13 



LCS includes the LAHET and HMCNP Monte-Carlo 
radiation transport codes. LAHET follows protons 
of all energies and secondary neutrons above 20 MeV 
throughout the geometry using a selection of on-line 
nuclear models and parameters; for these calculations 
LAHET2.83 was used with default selections except for 
use of ISABEL intranuclear-cascade model with de- 
fault parameters, prccquilibrium model following in- 
tranuclear cascade, and nuclear elastic scattering for 
neutrons and protons. 

LAHET passes the position, direction, and energy 
of all neutrons below 20 MeV as source particles to 
HMCNP, which uses evaluated cross-section data to fol- 
low all neutrons (and photons, if desired) throughout 
the geometry. 

The regional nuclide production (At. /proton) from 
LAHET and the binned neutron fluence (n/bin-p- 
cm 2 ) from HMCNP are scaled by region volume (cm 3 ) 
and proton source strength (p/s) to form production 
probabilities (At./cm 3 -s) and neutron fluxes (n/cm 2 -s). 
These are used with evaluated cross sections and decay 
data in CINDER'90 to calculate the temporal nuclide 
inventory during and following irradiation. 

In these simulations of thin-target measurements in 
environment having no features contributing to lower- 
energy neutron flux, contributions to nuclide inventory 
from the lower energy neutron reactions in HMCNP 
and in CINDER'90, although included, are negligible. 

To investigate the real predictive power of the codes 
but not just an ability to describe some known mea- 
sured data, all simulations have been performed at Los 



Alamos strictly before receiving any experimental data 
from Moscow that could be compared with theoretical 
results. 

PROCESSING EXPERIMENTAL AND SIMULATED 
DATA 

Determination of experimental activities and dose 
rates 

Activity of each nuclide in experimental samples 
was calculated via least square approximating nuclide 
counting rates of corresponding 7-cnergy peaks taking 
into account 7-yields and absolute effectiveness of the 
spectrometer at corresponding 7-energy. The total ac- 
tivities of the targets were defined as a sum of activ- 
ities of individual nuclides. The individual activities 
were determined dependently of properties of precur- 
sors present for practically all measured products: 

1) If precursors are absent or have lifetimes that are 
cither much less than the first measurement time after 
irradiation or much higher than the last measurement 
time, the following formula was used for activity calcu- 
lation: 

a (t,T) = — exp[-A(t - T)l . 
erj 

2) If precursors have lifetimes that are comparable 
with times passed after irradiation, the following for- 
mula was used: 

a{t,T) = — {A 1 exp[-X 1 (t-T)}+A 2 exp[-X 2 (t-T)}} , 

STj 

where, A , Ai, and A 2 are the coefficients determined 
via the least square approximation; A, Ai, and A 2 arc 
the decay constants of the nuclide and its precursors; 
t is the time after starting the irradiation; T is the 
irradiation duration; n is the 7-yield; e is the absolute 
effectiveness of the spectrometer. 

The dose rate at 1 cm distance from a target was de- 
fined via summing individual nuclides dose rates which 
were determined via individual activities multiplied by 
iVy-dose coefficients calculated for each nuclide using 
methods and data described in [5] . 

Simulated activities and dose rates 

To have a correct comparison of obtained experimen- 
tal results with the simulated ones, the simulated total 
activities were calculated by summing the individual 
simulated activities of only the nuclides that were iden- 
tified in measurements and used for calculating experi- 
mental activities. The simulated dose rates were calcu- 
lated by summing the individual dose rates of the same 
nuclides. 
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AGREEMENT BETWEEN SIMULATED AND EX- 
PERIMENTAL RESULTS 

The determined experimental and simulated activ- 
ities are presented in Figures 1-6. Additionally, the 
figures show the mean squared deviation factors calcu- 
lated via 

/ / / theor. value\\ 2 
(F) =loV V v exp ' value // 

for each of the data sets compared, where (...) desig- 
nates averaging over all of the comparison events. 

As seen from the figures, the simulated results are 
in good agreement with experimental data for most of 
the targets. In most cases, the mean squared deviation 
factor is from 1.15 to 1.47. The following variations of 
simulated-to-experimental agreements along time range 
are observed: 

1. nat Pb and 209 Bi: the simulated activities during 
the first 15 hours after irradiation are on the av- 
erage underestimated by a factor of 2, meanwhile, 
the results are in a much better agreement later. 

2. 63 Cu, E p = 1186 MeV: the simulated results, both 
for activity and dose rate, are in average 3 times 
underestimated during the first 5-300 hours after 
irradiation, meanwhile the discrepancies are negli- 
gible for other times. 

3. 63 Cu and 65 Cu, E p = 127 MeV: the simulated re- 
sults, both of activity and dose rate, are in aver- 
age 1.5 times lower for 63 Cu, and 2 times under- 
estimated for 65 Cu after ~ 1 day after irradiation, 
meanwhile during the first day after irradiation the 
results either practically coincide with experimen- 
tal data ( 63 Cu) or lie very close to them ( 65 Cu). 

The largest discrepancy between simulated and ex- 
perimental data is observed for 59 Co irradiated with 127 
MeV protons. Simulated activities are about 4 times, 
and dose rates - 2 times underestimated almost for the 
entire range of time. This reflects, probably, the im- 
portance of nuclear structure effects for this reaction 
(relatively low bombarding energy and vicinity of the 
target 59 Co to the doubly magic 56 Ni and magic 57 Ni 
nuclides) and an inadequate capability of LAHET to 
describe nuclide production near closed shells, which 
could be a result of using not good enough approx- 
imations for level densities, for inverse cross sections 
of evaporated and preequilibrium particles, and/or, for 
nuclear masses and binding energies near closed shells. 

CONCLUSIONS 

As a whole, the simulated activities and dose rates 
are in good agreement with experimental data for most 



of the targets along the whole time range. Neverthe- 
less, some serious discrepancies at certain times after 
irradiation, especially for 59 Co, were obtained. 

At the moment, we have not made a detailed investi- 
gation of the causes of the observed discrepancies. The 
list of products that are produced in the irradiations 
and/or are observed at decay phase and much con- 
tributes to the discrepancies observed have not been 
defined yet. Also, the contribution of various modes of 
nuclides formation in nuclear interactions (spallation, 
fission, fragmentation) to the observed discrepancies 
has not been analyzed properly. 

This work is our first step in the study of experi- 
mental and simulated activities of targets induced by 
intermediate and high energy protons. The analysis of 
observed discrepancies will be continued, and, if pos- 
sible, will be extended to thick targets in which sec- 
ondary neutrons, rather than the incident protons, have 
the largest contribution to target activation. The use 
of the ISABEL model in LAHET for the higher-energy 
simulations above its nominal 1-GeV limit could lead to 
inconsistencies; similar calculations using the alternate 
Bertini model, or preferably the nuclear models of the 
CEM95 code [6], in LAHET are called for. 

The experimental part of this work have been sup- 
ported by the ISTC Project #839 and the theoretical 
study was supported by the U. S. Department of Energy 
under contract W-7405-Eng-36. 
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Figure 1: Experimental activities and dose rates and results calculated with the LCS-CINDER code system for nat Pb 
and 209 Bi irradiated with 1478 MeV protons. The filled symbols show calculated values during the irradiation. The 
mean squared deviation factors (see text) are shown for both. 
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Figure 2: The same as Fig. 1, but for 59 Co irradiated with 1186 MeV and 127 MeV protons. 
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Figure 3: The same as Fig. 1, but for 63 Co and 65 Co irradiated with 1186 MeV protons. 
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Figure 4: The same as Fig. 1, but for 63 Co and 65 Co irradiated with 127 MeV protons. 
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Figure 5: The same as Fig. 1, but for 232 Th irradiated with 97 MeV and 795 MeV protons. 
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Figure 6: The same as Fig. 1, but for nat U irradiated with 97 MeV and 795 MeV protons. 
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